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Chemical  Reaction  in  a  Glassy  Matrix: 

Dynamics  of  Ligand  Binding  to  Protohcmc 
in  Glycerol:  Water 

J.R.  Hitt1,  M.J.  Cote1,  C  D.  Dlott1,  J.F.  Kiuffmsn' ,  J.D.  McDonald* . 
P.J.  Steinbtch J.R.  Berendien1 ,  and  H.  Fnuenfelder * 

'School  of  Chemical  Sciences,  University  of  Illinois  at  Urbane-Champaign, 
Urbana,  1L  61801,  USA 

department  of  Physics,  University  of  Illinois  at  Urhana-Champaign, 
Urbana,  1L  61801,  USA 


A  fundamental  understanding  of  condensed  phase  chemical  reaction  dynamics 
can  be  obtained  from  the  study  of  biomolecules,  particularly  h«e  and  heme- 
proteins.  The  complexity  of  these  systems  gives  rise  to  a  rich  variety  of 
phenanena,  allowing  many  aspects  of  condensed  Batter  reactions  to  be 
examined.  The  rate  theories  and  puzzles  of  hemeproteln  kinetics  have 
recently  been  discussed  by  Frauenfelder  and  Wolynes  [  1 ).  In  this  work  we 
present  new  experiments  on  ligand  binding  to  protoheme  (Fe: protoporphyrin- 
IX)  in  a  glassy  matrix.  We  have  studied  the  rebinding  of  carbon  monoxide 
arter  photodissociation  over  a  wide  range  of  time  [5ps-l0ss]  and  temperature 
[300K-70K].  The  significance  of  our  results  is  that  (1)  the  influence  or 
friction  and  nonadlbatlclty  on  a  condensed  phase  reaction  can  be  directly 
investigated.  (2)  the  Influence  of  an  Inhomogeneous  glassy  matrix  (glycerol* 
water  75:25)  on  the  reaction  can  be  studied,  and  (3)  meaningful  comparison 
between  protoheme  and  heneproteln  kinetics  Isolates  the  role  or  the  protein 
relaxation  in  the  reaction. 

Tmsperature- dependent  ligand  rebinding  to  protoheme  was  studied  with  us 
resolution  by  Alberdlng  et  al.  (2).  We  extended  the  results  to  the  pa  time- 
scale  using  conventional  pump-probe  methods  with  565  nm  dye  laser  pulses 
OOpJ,  3ps,  FWHH.  80Hz).  Figure  1  shows  a  log-log  plot  of  the  time- resolved 
absorbance  change  (proportional  to  H(t),  the  fraction  of  protoheme  molecules 
that  have  not  rebound  CO  at  time  t)  at  300K  and  260K.  The  slow  process 
(Process  S)  Involves  binding  of  CO  that  have  escaped  their  initial  solvent 
cage  by  diffusion.  This  solvent  process  is  exponential  In  lime  with  a 
pseudo-first-order  rate  coefficient  that  depends  on  [CO).  As  T  Is  de¬ 
creased,  the  fraction  or  CO  molecules  that  escape  the  cage  decreases  until 
it  is  below  0.01  at  2*0K  (2),  well  above  the  glass  transition  Tc  •  173*. 
Subtraction  of  the  260K  Process  S  data  (solid  circles  from  [2))  yields  a 
nonexponential  decay  (Process  1),  the  solid  rectangles  in  Fig.  1.  The 
relative  fraction  of  ligands  that  rebind  via  Process  1  Increases  as  T  de¬ 
creases  from  300K  to  2*ok.  Process  I  is  nonexponential  at  all  temperatures 
up  to  300K. 

As  T  Is  decreased  below  2*0K,  process  S  vanishes,  the  rate  of  process  1 
decreases  in  a  manner  consistent  with  [2),  and  a  new  process.  1*.  Is 
observed.  The  amplitude  of  1*  Increases  with  farther  decrease  In  T  until  it 
dominates  1.  Figure  2  shows  the  rebinding  data  at  U0K  (open  rectangles), 
the  temperature  at  which  processes  Z  and  X*  have  equal  amplitudes.  Subtrac¬ 
tion  of  process  1,  which  is  kn**n  to  have  the  form  of  a  power-law  decay, 
reveals  that  process  Ic  (open  circles)  decays  exponentially  with  a  15  ps 
lifetime.  As  T  is  decreased  from  220  to  70X,  the  rate  of  process  X*  da- 
creases  linearly  with  T. 
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Fig.  2.  Ligand  rebinding  data  at  1*OK  where  the  fraction  of  Uganda  that 
rebind  via  process  1  la  equal  to  the  fraction  that  rebind  via  process  I*. 


Figure  3  svawnarlxes  the  twperature  dependence  of  ligand  rebinding  via 
these  three  processes.  At  the  lowest  temperatures  studied,  process  I* 
dominates.  As  the  taeperatire  increases,  the  fraction  of  ligands  rebinding 
via  process  1  Increases  until  2«0K.  At  that  point  the  solvent  process 
Increases  In  Importance. 

In  [2],  process  1  was  interpreted  to  arise  from  recombination  of  Uganda 
trapped  within  the  solvent  cage.  The  barrier  for  this  process  Is  a  property 
of  the  hene,  and  probably  arises  as  a  consequence  of  the  transition  frem  the 
unbound  S-2  domed  heee  structure  to  the  bound  $-0  planar  structure.  The 
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rebinding  is  nonexponent i«i  because  this  barrier  possoncs  distributed 
activation  enthalpy.  In  the  6 -chain  of  hemoglobin,  process  1  becomes  nearly 
exponential  above  ca.  250K,  because  the  protein  structure  relaxes 
significantly  on  the  timescale  (t  ns)  of  tne  rebinding  lr.  prptones-e, 

this  conformational  relaxation  does  not  occur  on  the  fester  liptscale  t*»0 
ps)  of  the  rebinding. 

The  observed  (linear)  temperature  dependence  and  the  large  rate  of 
process  1*  at  lot#  temperatures  suggest  that  it  does  not  occur  ty  classical 
over' the- bar rler  motion.  The  exponential  time  dependence  shows  that  the 
height  of  tne  barrier  is  not  distributed.  Process  t*  e*y  involve  recam* 
blnatlon  to  sn  S-1  heme  which  has  not  dosed,  At  t*»0K,  the  rebinding  rale  of 
TEtO/s  is  a  factor  of  -*0  slower  than  *T/n,  the  frequency  or  collision 
between  the  ligand  and  the  heme.  The  reduction  of  the  recombination  rate 
can  be  caused  by  entropy,  friction,  and  nonadiabatic  barrier  crossing  fl). 

This  research  was  supported  In  part  by  the  National  Science  foundation. 
Solid  State  Chemistry,  grant  DMft  8A-15070  <000),  the  U.  S.  Qe;.jrtstr,t  of 
Health  and  Human  Services  under  Grant  Gtt  18051,  and  the  National  Science 
Foundation  under  Grant  Df©  82-09616  (HF),  and  the  National  Science  Founds* 
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INTRODUCTION 

Protein  motions  are  studied  using  the  binding  of  carbon  monoxide  (CO)  to  myo¬ 
globin  (Mb)  as  probe.  A  number  of  experimental  tools  have  been  used  to  study  this 
simple  protein  reaction  (1,2),  but  the  most  detailed  information  has  come  from 
flash  photolysis  experiments  over  wide  ranges  in  time,  temperature,  viscosity,  and 
pressure  (3,4). 

Here  I  describe  the  results  of  flash  photolysis  experiments  using  “pressure 
titration"  to  probe  large-scale  protein  motions  (5,6).  In  a  pressure  titration 
experiment,  an  ensemble  of  proteins  is  prepared  in  different  nonequilibrium  states 
by  using  different  paths  in  the  pressure-temperature  plane  to  reach  a  given  pres¬ 
sure  and  (low)  temperature.  In  the  case  of  CO  binding  to  Mb  the  relaxation  of  the 
recombination  kinetics  between  two  nonequilibrium  states  of  the  protein  ensemble 
at  the  same  temperature  and  pressure  probes  large-scale  protein  motions.  These 
motions  have  been  previously  denoted  as  EF1  (equilibrium  fluctuations)  and  FIM1 
(^functionally  Jmportant  motions)  (5,7).  The  pressure  titration  experiments  imply 
that  the  EF1  must  be  described  by  a  distribution  of  activation  energies  ranging 
between  70  to  85  kj/mol. 


HIERARCHICAL  PROTEIN  MODEL 


The  principal  features  of  the  hierarchical  model  of  proteins  are  (7):  (1) 

A  globular  protein  exists  in  a  large  number  of  conformational  substates  (CS) 
separated  by  barriers,  some  of  which  become  effectively  infinitely  high  below 
about  160  K  (CS1).  (11)  Transitions  among  the  substates  are  equilibrium  fluctua¬ 

tions  (EF);  nonequilibrium  transitions  from  one  protein  state,  say  MbCO,  to 
another  state,  say  deoxyMb,  Involve  functionally  Important  motions  (FIMs).  (Ill) 
Substates  and  equilibrium  fluctuations  possess  hierarchies.  In  Mb  four  different 
tiers  of  substates  and  motions  are  observed,  (lv)  The  different  tiers  of  sub¬ 
states  are  described -by  distributions  of  relaxation  rates. 


The  pressure  titration  experiments  probe  equilibrium  fluctuations  (EF1)  in 
tier  1  of  the  hierarchy  of  conformational  substates  (CS1).  The  EF1  most  likely 
involve  large-scale  motions  In  the  protein  and  possibly  Include  the  hydration 
layer  and  part  of  the  solvent.  FIM1  may  be  the  motion  that  permits  ligands  to 
enter  and  leave  the  protein. 


PRESSURE  TITRATION  EXPERIMENTS 


Below  about  180  Kt  Mb  rebinds  CO  after  photodissociation,  but  the  kinetics 
follow  a  power  law  in  time.  This  means  that  the  protein  cannot  be  in  a  unique 
state  but  must  exist  in  a  large  number  of  conformational  substates  (CS1)  with 
different  barriers  for  recombination  from  the  heme  pocket.  The  fraction  N(t)  of 
proteins  which  have  not  rebound  a  ligand  at  time  t  after  photodissociation  is 
given  by 

N(t)  -  / g(H)  exp(-k(H)t)dH  (1) 

where  g(H)dH  is  the  probability  of  a  protein  having  a  barrier  with  activation 

•  » 

enthalpy  between  H  and  H  +  dH  and  k(H)  is  the  rate  coefficient.  If  rebinding 
is  thermally  activated  then  k(H)  ■  Aexp(-H/RT).  The  distribution  g(H)  and  pre¬ 
exponential  A  are  found  from  the  experimentally  measured  N(t).  Alternatively  I 
use  the  model  of  Young  and  Bowne  for  the  distribution  g(H)  which  yields  excel¬ 
lent  fits  to  the  low  temperature  kinetics  (A, 8). 

Under  pressure  the  binding  of  CO  to  Mb  speeds  up  (6,9).  The  effect  is 
particularly  dramatic  and  relatively  simple  to  interpret  at  low  temperatures. 

The  ensemble  of  proteins  is  prepared  in  different  nonequilibrium  states  as 
follows:  High  pressure  (2  kbar)  is  applied  to  the  sample  at  high  temperature 
(300  K)  and  the  sample  is  then  cooled  to  say  100  K  before  the  pressure  is  re¬ 
leased  to  1  bar  (state  6').  The  sample  can  also  be  cooled  to  low  temperature 
(100  K)  at  atmospheric  pressure  (state  8).  Since  CO  binding  is  much  faster  at 
state  S'  than  at  state  B,  the  two  states  are  not  in  equilibrium  and  can  be 
distinguished  by  measuring  N(t).  The  g(K)  for  state  B'  has  a  similar  shape  as 
g(H)  for  state  8  but  is  simply  shifted  to  smaller  actlvstion  enthalpies. 

The  difference  in  the  distributions  g(H)  is  explained  simply  in  terms  of 
the  conformational  substates  of  tier  1.  At  6  the  ensemb'e  of  proteins  is  in  a 
distribution  of  CS1  characteristic  of  1  bar  while  at  B'  the  distribution  of  CS1 
is  characteristic  of  2  kbar.  Transitions  between  states  B'  and  B  correspond  to 
EF1.  The  EF1  are  induced  as  follows:  Starting  in  state  B'  at  say  100  K  and 
1  bar,  the  recombination  kinetics  N(t)  are  measured.  The  temperature  is  then 


Increased  to  Tm,  left  there  for  600s  (the  annealing  time,  ta),  cooled  again  to 
100  K  and  N(t)  measured.  This  procedure  Is  repeated  with  increasing  Ta  until 
N(t)  has  relaxed  to  follow  the  curve  corresponding  to  state  0.  The  result  Is 
definitive:  Op  to  Ta  *  160  K,  N(t)  for  6*  is  unrelaxed  but  at  Ta  »  200  K  N(t) 
has  completely  relaxed  Co  the  kinetics  corresponding  to  0. 

Since  the  distributions  g(H)  for  0  and  0'  have  the  same  shape,  the  states 
0  and  0'  are  characterized  by  the  minimum  activation  enthalpy  Hm^n  of  the  dis¬ 
tributions  (0).  The  extent  of  relaxation  is  characterized  by  the  fractional 
shift  0(ta,To)  of  from  1^(0):  0(ta,Ta)  -  l^ln(t..T11)-HBln(0) J/ 

(Hain(0' )-UBin(0) ] .  Assuming  that  the  relaxation  0*  ♦  0  is  described  by  a 
distribution  of  activation  energies,  A(ta,Tm)  is  given  by 

6(ta,To)  -  /p(Ej)  exp(-kjta)dE,  (2) 

where  p(E})dE1  is  the  probability  of  finding  a  barrier  between  Ej  and  Ej  +  dEp 
l  fit  the  data  with  kj  ■  A^exp(-Ej/RT)  and  a  box  distribution, 

E  -E  ,  )”1,  E  ,  <  E,  <  E 

max  min  sun  1  max 

•  (3) 

0  ,  otherwise 

•« 

A  good  fit  to  the  relaxation  of  N(t)  parameterized  by  Hmin  is  obtained  for 
Aj  “  ID1*'3  a"1,  E^n  “  70  kJ/mol,  and  Emax  -  85  kJ/mol. 

DISCUSSION 

The  pressure  titration  experiments  on  CO  rebinding  to  Mb  after  photodisso¬ 
ciation  show  that  at  1  bar  the  EF1  set  in  at  160  K  and  extend  over  a  range  of 
about*  A0  K.  The  resulting  distribution  of  relaxation  rates  is  reminiscent  of 
glassy  relaxation  processes  (10,11,12).  It  is  possible  that  the  hydration  layer 
and  solvent  strongly  influence  the  motions  of  tier  1.  Additional  evidence  for 
EF1  has  also  come  from  Raleigh  scattering  of  MSssbauer  radiation  (13)  and 
microwave  absorption  in  the  hydration  shell  (1A).  However  we  are  only  at  the 
beginning  of  the  study  of  EF1  and  FIM1.  Clearly  more  theory  is  needed  and  more 
experiments  are  required  to  explore  the  properties  of  EF1  and  FIM1 ,  to  study  the 
role  of  the  hydration  layer  and  solvent  (15),  and  to  determine  the  connection 
between  EF1  and  unfolding. 
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iron.  The  CO  molecule  executes  a  Brownian  motion  in  the  solvent,  motions  of  the  protein  into  account.  In  a  second  model  we  use 

moves  into  the  protein  matrix  M,  migrates  through  the  matrix  into  ligand  binding  as  a  perturbation  that  induces  protein  motions  and 

the  heme  pocket  B  and  finally  binds  covalently  to  the  heme  iron.  look  at  these  motions  alone.  The  perturbation  caused  by  ligation 

In  the  one-particle  model, 12  we  assume  that  the  protein  forms  or  delegation  leads  to  a  chain  of  events  with  many  similarities  to 
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discussed  CO  binding  and  concentrate  on  the  proteinquake  here. 
dHeAg(HBA)  exp(  -  Hb  / HT)  (7)  Experimental  observations  to  be  discussed  below  suggest  that 

S'A/  '  BA  the  relaxation  of  MbCO  after  photodissociation  occurs  through  at 


of  protein  substates  given  in  Fig.  12. 25  A  given  state,  for  instance 
MbCO,  shown  at  the  top  of  Fig.  12  as  a  single  potential  well,  is 
subdivided  into  a  large  number  of  conformational  substates  (CS')< 
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ABSTRACT 

The  Infrared  stretching  bands  of  carboxymyoglobin  (MbCO)  and  the 
rebinding  of  CO  to  Mb  after  photodissociation  have  been  studied  in  the 
temperature  range  from  10  to  300K  in  a  variety  of  solvents.  Four  stretching 
bands  imply  that  MbCO  can  exist  in  four  substates,  Aq  to  Aj.  The  temperature 
dependences  of  the  intensities  of  the  four  bands  yield  the  relative  binding 
enthalpies  and  entropies.  The  integrated  absorbances  and  the  pH  dependences 
of  the  bands  permit  identification  of  the  substates  with  the  conformations 
observed  in  the  X-ray  data  [Kuriyan  et  al. ,  J.  Mol.  Biol.  192  (1986)  133).  At 
low  pH,  Aq  is  hydrogen-bonded  to  His  E7.  The  substates  Aq  to  Aj  interconvert 
above  about  180K  in  a  752  glycerol-water  solvent  and  above  270K  in  buffered 
water.  No  major  Interconversion  is  seen  at  any  temperature  if  MbCO  is  embedded 
in  a  solid  poly (vinyl )alcohol  matrix.  The  dependence  of  the  transition  on 
solvent  characteristics  is  explained  as  a  slaved  glass  transition. 

After  photodissociation  at  low  temperature  the  CO  is  in  the  heme  pocket 
B.  The  resulting  CO  stretching  bands  which  are  identified  as  B  substates  are 
blue-shifted  from  those  of  the  A  substates.  At  40K,  rebinding  after  flash 
photolysis  has  been  studied  in  the  Soret,  the  near-XR,  snd  the  integrated  A 
and  B  substates.  All  data  lie  on  the  same  rebinding  curve  and  demonstrate 
that  rebinding  is  nonexponential  in  time  from  at  least  100  ns  to  100  ks.  No 
evidence  for  discrete  exponentials  is  found. 

Flash  photolysis  with  monitoring  in  the  IR  shows  four  different  pathways 
within  the  pocket  B  to  the  bound  substates  A^ .  Rebinding  in  each  of  the  four 
pathways  B  ♦  A  is  nonexponential  in  rime  to  at  least  10  ks  and  the  four 
pathways  have  different  kinetics  below  180K.  From  the  time  and  temperature 
dependence  of  the  rebinding,  activation  enthalpy  distributions  g(HBA)  and 
preexponentials  *BA  are  extracted.  No  pumping  from  one  A  substate  to  another, 
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or  one  B  substate  to  another,  is  observed  below  the  transition  temperature  of 
about  180K. 

If  MbCO  is  exposed  to  intense  white  light  for  10  to  10^s  before  being 
fully  photolyzed  by  a  laser  flash,  the  amplitude  of  the  long-lived  stares 
increases.  The  effect  is  explained  in  terms  of  a  hierarchy  of  substates  and 
substate  symmetry  breaking. 

The  characteristics  of  the  CO  stretching  bands  and  of  the  rebinding 
processes  in  the  heme  pocket  depend  strongly  on  the  external  parameters  of 
solvent,  pH,  and  pressure.  This  sensitivity  suggests  possible  control 
mechanisms  for  protein  reactions. 
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INTRODUCTION 


The  binding  of  smell  ligands  such  as  carbon  monoxide  (CO)  to  heme 
proteins  at  first  appears  to  be  a  simple  one-step  process1,  but  a  closer  look 
over  extended  ranges  in  temperature,  time,  and  wavelength  reveals  a  surprising 

complexity.  Binding  is  governed  by  a  sequence  of  barriers  rather  than  a 
2 

single  one.  After  photodissociation,  rebinding  occurs  at  temperatures  as  low 
3 

as  a  few  K  and  is  nonexponential  in  time  below  180K  in  a  752  glycerol-water 
solvent .2’“* ’ "*  Below  about  40K,  binding  no  longer  proceeds  by  classical  over- 
the-barrler  motion,  but  through  quantum-mechanical  tunneling. The  behavior 

of  the  photodlssociated  CO  within  the  heme  pocket,  followed  by  IR  spectros- 

8  9 

copy,  turns  out  to  be  complicated.  ’ 

In  the  present  paper  we  investigate  processes  within  the  heme  pocket  of 
myoglobin  (Mb)  in  more  detail.  Mb  has  been  called  the  hydrogen  atom  of 
biology;  it  is  a  globular  protein  of  about  18  kdalton  that  stores  dioxygen 
(02)  in  muscles.10  02  or  CO  binds  reversibly  to  the  central  iron  atom  of  the 
heme  group  in  Mb.  Since  Mb  is  sturdy  and  very  well  explored,  it  offers  a 
nearly  ideal  system  for  more  sophisticated  studies. 

The  impetus  for  the  present  work  came  from  two  directions:  (1)  At  low 

8  9 

temperatures,  CO  bound  to  Mb  shows  a  number  of  stretching  bands.  ’  During 
our  earlier  work  we  had  noticed  that  the  different  bands  rebind  with  different 
rates  after  photodissociation.  A  rate  difference  has  also  been  noticed  by 
Gerwert  et  al.  and  M.  Chance  et  al.  Such  a  rate  variation  suggests  a 
control  mechanism  and  calls  for  more  studies.  (2)  Nonexponential  kinetics 
below  180K  has  been  seen  in  all  protein-ligand  systems  that  have  been  studied. 
In  general,  nonexponential  time  dependences  can  be  explained  either  by 
homogeneous  or  Inhomogeneous  processes.  By  using  a  sequence  of  laser  flashes 
repeated  before  all  tie  Mb  molecules  have  rebound  a  CO  ("hole  burning  in 
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Cine”)  we  have  demonstrated  that  puaplng  f ron  short-lived  to  long-lived  states 
at  60  and  70K  is  very  saall  and  that  the  nonexponential  kinetics  results 
predoalnantly  from  lnhoaogeneitles  among  Mb  molecules  and  cannot  be  explained 
by  multiple  states  within  one  protein  molecule.^'1^’1*  In  view  of  these 
results  it  is  surprising  that  B.  Chance  et  al.  recently  claimed  evidence  for  a 
bl-  or  triexponential  time  course  and  considerable  puaplng  froa  short-  to 
long-lived  components  in  the  binding  of  CO  to  Mb  at  40K.^  The  difference 
between  their  interpretation  and  ours  is  fundamental.  B.  Chance  et  al.  assume 
that  the  different  lifetimes  are  caused  by  different  pathways  or  sites  within 
the  same  molecule.  We  postulate  a  large  number  of  structurally  different 
conformational  substates  (CS),  with  each  molecule  remaining  In  a  given  CS  at 
temperatures  where  the  solvent  is  frozen.^’16  Our  Interpretation  leads  to  the 

concept  of  a  hierarchy  of  substates17  and  to  connections  between  proteins  and 
1 8 

(spin)  glasses.  We  will  discuss  the  concept  of  CS  and  of  the  hierarchy  of 
CS  in  Subsection  4.5. 

19 

Here  we  report  on  expanded  low-temperature  experiments.  The  new 
experiments  verify  our  earlier  experiments  and  also  yield  surprising  results 
that  open  new  windows  on  phenomena  in  the  heme  pocket. 
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2.  EXPERIMENTAL  APPROACH 

The  primary  experimental  technique  la  low-temperature  flash  photolysis. 

An  MbCO  sample  at  a  controlled  temperature  Is  photodlssoclated  with  a  light 
pulse  and  rebinding  is  monitored  through  optical  absorption  measurements.  We 
have  observed  rebinding  at  many  different  wavelengths,  as  Indicated  In  Fig.  1, 

thus  monitoring  the  rebinding  process(es)  at  different  locations  In  the  heme 
20  21 

pocket.  ’  In  the  Soret  and  the  a  and  6  bands,  changes  in  the  electronic 
structure  of  the  heme  on  binding  are  seen.  The  near-IR  band  (band  III,  near 
760  nm)  is  a  charge-transfer  band  and  Involves  both  the  Iron  atom  and  the 

heme.  In  the  IR,  the  bands  near  1950  cm  1  result  from  the  stretching  mode  of 
22  23  “1 

the  bound  CO  ’  and  the  ones  near  2100  cm  result  from  the  stretching  mode 
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of  the  photodlssoclated  CO  within  the  pocket.  ’  If  only  one  simple  rebinding 
process  were  to  occur,  all  of  these  markers  would  give  the  same  time  course. 

Four  different  systems  were  used  for  the  experiments,  a  nanosecond  (ns) 
visible-wavelength  flash  photolysis  system,  a  conventional  dispersive  vislble- 
near-IR  spectrophotometer,  an  FTIR  spectrophotometer  and  a  microsecond  mid-IR 
flash  photolysis  system.  Samples  were  prepared  from  lyophilized  sperm  whale 
Mb  from  Sigma  (St.  Louis,  M0)  dissolved  in  buffer  solution.  Buffers  used 
were:  0.4M  citrate-phosphate  buffer  for  the  pH  5.0  sample,  0.4M  carbonate- 
bicarbonate  buffer  for  pH  9.0  samples,  0.1M  HEPES  buffer  for  the  pH  7.0 
poly(vinyl)alcohol  (PVA)  samples,  and  0.5M  sodium  phosphate  buffer  for  all  of 
the  other  samples.  Solvents  were  water,  75Z  (volume /volume)  glycerol-water, 
and  PVA.  The  Soret  data  were  taken  on  10-14  uM  samples  placed  in  a  1  cm  x  1 
cm  cuvette  and  the  near-IR  experiments  were  done  on  a  1.5  mM  sample  in  a  1  cm 
by  0.1  cm  cell. FTIR  and  fast-IR  measurements  were  performed  on  15-20  mM 
samples  held  between  two  CaFj  windows  with  a  75  urn  mylar  spacer. 
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The  ns  system  employed  a  30  ns  pulse  from  a  frequency-doubled,  Q-switched 
Nd+-glass  laser  (530  nm,  300  mJ).  Rebinding  was  monitored  with  the  light  from 
a  tungsten  lamp  passed  through  a  monochromator.  Monitor  beam  Intensities  were 
kept  small  (<  10  photons/s  absorbed  per  protein  molecule).  The  photo¬ 
multiplier  (Hamamatsu  928)  signal  was  digitized  with  a  LeCroy  TR8818  transient 
digitizer  from  10  ns  to  300  ys  and  with  a  logarithmic  time-base  digitizer  from 
2  ps  to  300  s.  A  storage  cryostat  (Janis  10-DT)  permitted  temperature 
control. 

Our  OLIS-Cary  14  spectrophotometer  is  interfaced  to  an  IBM  PC/AT  computer. 
Kinetics  data  obtained  from  this  instrument  extend  from  60  to  3  *  10^s. 

The  FTIR  experiments  were  performed  on  a  Mattson  Sirius  100  FTIR  spec¬ 
trometer.  All  spectra  were  taken  at  1  or  2  cm  *  resolution.  Low  temperatures 
were  produced  with  a  closed-cycle  helium  refrigerator,  CTI  model  21. 

Photolysis  in  the  kinetics  experiments  was  achieved  with  a  500  ns,  0.3  J, 

590  nm  pulse  from  a  Phasu-R  DL2100C  dye  laser  using  rhodamlne  6G  dye. 

Absolute  absorbance  spectra  were  obtained  by  subtracting  the  solvent  back¬ 
ground  spectrum  from  the  MbCO  spectrum.  The  difference  absorption  spectra  in 
the  kinetics  experiments  were  obtained  by  referencing  the  photolyzed  spectrum 
to  the  unphotolyzed  background  spectrum. 

The  fast  kinetics  of  the  CO  stretching  bands  was  measured  in  a  separate 
mid-IR  flash  photolysis  system.  The  Infrared  monitoring  light  was  produced  by 
a  Spectra-Physics/Laser  Analytics  tunable-diode  laser  and  directed  through  the 
protein  solution  and  monochromator  onto  a  U^-cooled  IR-Associates  HgCdTe 
detector.  Photolysis  was  induced  by  the  same  laser  usad  in  the  FTIR  experi¬ 
ments.  For  times  between  100  ms  and  100  s  the  diode  laser  current  was  modu¬ 
lated  at  10  kHz  and  the  signal  was  amplified  by  a  PAR  5101  lock-in  amplifier. 
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3.  RESULTS  AND  DATA  EVALUATION 

A  coherent  picture  of  the  low-temperature  binding  emerges  after  the  data 
obtained  with  the  different  Instruments  and  in  the  different  wavelength 
regions  are  pieced  together. 

3.1  IR  Measurements  of  MbCO 

(I)  The  A  Substates.  The  stretching  bands  of  the  bound  CO  are  superb 
probes.  He  have  measured  the  "dark"  or  bound-state  IR  spectrum  of  MbCO  at 
temperatures  between  15  and  300K  in  75T  glycerol-water  (Fig.  2),  in  water 
(Fig.  3),  and  in  solid  poly(vinyl)alcohol  (PVA)  (Fig.  4). 

The  spectra  shown  in  Figs.  2-4  can  be  decomposed  into  Gaussian  super- 

24  25 

positions  of  Lorentzians,  called  Volgtlans.  ’  Two  such  decompositions  are 
shown  in  Fig.  5.  Four  different  IR  bands  can  be  distinguished.  We  denote 
the  bands  in  order  of  decreasing  wavenumbers  by  Aq  ("  1966  cm-1),  Aj  (“  1946 
cm-1),  A2  (*  1941  cm-1),  and  A-j  («  1930  cm-1).  (Because  of  the  new  band  Aj 

Q 

the  nomenclature  differs  from  the  one  we  used  earlier.  )  The  fits  yield  the 
center  positions  v^q,  the  full-width  at  half  maximum  T,  and  the  areas  of  the 
four  bands.  Since  the  band  ^  1®  difficult  to  extract  unambiguously,  we  shall 
discuss  only  the  other  three. 

(II)  Temperature  Dependence  of  the  A  Bands.  In  order  to  characterize  the 
temperature  dependence  of  the  individual  bands  in  glycerol-water,  we  select 
the  most  intense  band,  Aj ,  as  standard  and  plot  in  Fig.  6  the  logarithm  of  the 
intensity  ratios  Aq/Aj  and  Aj/Aj  as  a  function  of  103/T.  The  ratios  A}/Aj  in 
Fig.  6  show  a  different  behavior  for  PVA  and  glycerol-water.  For  PVA,  the 
ratios  change  very  little.  MbCO  in  75Z  glycerol-water,  in  contrast,  shows  two 
regions.  Above  about  180K,  the  ratios  change  rapidly;  oelow  that  temperature, 
the  change  is  much  slower.  The  behavior  above  180K  can  be  understood  if  the 
four  bound  substates  Aq  to  Aj  have  different  binding  enthalpies  and  entropies 
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and  may  interconvert  freely.  Below  180K,  the  change  may  be  explained  by 
assuming  chat  the  proteins  are  frozen  into  particular  substatet  and  the 
extinction  coefficients  for  the  different  A  substates  depend  slightly 
differently  on  temperature.  Above  180K,  the  ratios  are  given  by 

Ai/A1  -  exptSj/R]  expl-Hj/RTKCjCD/Ejd))  (1) 

wtv  and  Hj  are  the  entropy  and  enthalpy  of  substate  A^  with  respect  to 

ate  Aj  and  £j(T)  is  the  extinction  coefficient  of  substate  Aj  at 
temperature  T.  Assuming  as  a  first  approximation  that  the  ratio  of  extinction 
coefficients  is  unity,  relative  binding  enthalpies  and  entropies  have  been 
extracted  from  the  data  in  Fig.  6.  The  binding  parameters  for  two  values  of 
pH  are  given  in  Table  1. 

The  total  integrated  area  ZAj  of  the  A  bands  of  MbCO  in  751  glycerol- 
water,  given  as  a  function  of  temperature  in  Fig.  7,  shows  a  decrease  of  about 
20X  in  going  from  10  to  300K.  Such  a  decrease  can  be  explained  by  a  change  in 

It 

the  electron-nuclear  coupling. 

The  temperature  dependence  of  the  center  frequency  v^q  is  given  in  Fig.  8 
for  the  bands  A^  and  Aj  for  two  samples.  The  temperature  dependence  of  the 
line  width  T  is  displayed  in  Fig.  9. 

(lit)  Nonequi  li  brium  Behavior  of  the  A  Bands.  The  ratios  Aq/Aj  and  A^/Aj 
in  Fig.  6  follow  a  van't  Hoff  behavior  from  260K  to  190K.  Below  190K,  they 
deviate  from  the  van't  Hoff  line  and  become  essentially  temperature  inde¬ 
pendent.  We  explain  this  temperature  dependence  by  assuming  that  the  A  sub¬ 
states  are  not  in  dynamic  equilibrium  below  about  180K  in  75X  glycerol-water. 
To  verify  this  assumption  we  have  measured  the  time  dependence  of  the  IR 
spectrum  after  the  temperature  was  lowered  from  195K  to  the  measurement 
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temperature  T.  We  denote  with  Aj(t)  the  integrated  absorbance  of  band  Aj  at 
the  time  t  after  lowering  the  temperature.  Aj(0)  is  the  integrated  absorbance 
at  19SK,  and  A^(»)  the  equilibrium  value  at  temperature  T  which  would  be 
reached  after  an  extremely  long  («•)  time.  We  eatimate  Aj(~)  by  extrapolating 
the  van't  Hoff  line,  obtained  from  260  to  195K,  to  the  temperature  T.  The 
fraction  nj(t)  of  MbCO  molecules  that  have  not  yet  transferred  to  substate  Aj 
at  time  t  after  the  change  in  T  is  given  by 

n,(t)  -  {A^O-AjC-imAjW-A^-)}.  (2) 

ng(t)  at  190K  is  shown  in  Fig-  10.  The  data  Indicate  that  the  approach  to 
equilibrium  is  nonexponenr lal  in  time.  In  addition  we  observe  that  the  time 
to  reach  equilibrium  is  extremely  long  below  180K.  The  A  substates  are  not  in 
equilibrium  below  180K. 

22 

(lv)  External  Influence  on  A  Substates.  Makinen  et  al.  and  Brown 
27 

et  al.  have  pointed  out  that  external  Influences  can  change  the  relative 
intensities  of  the  A  substates.  Such  effects  can  be  seen  in  Figs.  2-4.  Aq  is 
much  more  intense  in  PVA  and  in  ice  than  in  glycerol-water.  Pressure  and  pH 
also  affect  the  ratio  Aj /Aj .  Fig.  11  shows  that  A^  decreases  while  Aj  and  Aq 
increase  with  increasing  pressure.  Fig.  12  shows  that  Aq  increases  at  the 
expense  of  Aj  with  decreasing  pH.  While  the  band  intensities  are  affected 
strongly  by  solvent  state,  pressure,  and  pH,  the  band  positions  remain  nearly 
unchanged. 

The  close  coupling  of  the  environment  with  the  protein  is  supported  by 
the  data  in  Fig.  13.  The  IR  spectrum  of  a  75X  glycerol-water  solvent  exhibits 
a  transition  near  180K,  essentially  at  the  temperature  where  the  substate 


exchange  is  frozen 
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3.2  Overall  Rebinding  at  40K 

Figure  1  shows  Chat  rebinding  afcer  photodissociation  can  be  monitored  at 
various  locations  In  the  heme-CO  system.  We  have  studied  the  rebinding  at  40K 
over  a  wide  range  of  does  at  various  wavelengths.  We  selected  40K  because  B. 
Chance  et  al.*"’  have  claimed  evidence  for  discrete  exponentials  at  this 
temperature.  The  results  of  our  measurements  are  shown  In  Fig.  14.  The  data 
correspond  to  four  different  spectral  regions.  The  symbol  Q  indicates 
measurement  at  440  nm,  near  the  deoxy  Soret  peak;  Q  indicates  the  Integrated 
absorbance  of  the  "760  nm"  charge  transfer  band  III;  V  refers  to  the  integral 
over  all  A  (bound)  CO  Btretching  bands;  and  the  0  shows  the  Integrated  magni¬ 
tude  of  the  absorbance  change  of  the  Soret  band  region  from  400  to  450  nm. 

The  data  In  Fig.  14  are  smooth  over  more  than  twelve  orders  of  magnitude 
In  time.  No  evidence  for  the  bi-  or  t ri exponent i al  behavior  is  visible. 

Within  experimental  uncertainties  all  of  the  above  measurements  give  the  same 
time  course. 

3.3  Rebinding  to  the  A  Substates 

New  features  appear  when  the  rebinding  of  each  of  the  four  A  substates  is 
monitored  separately.  Fig.  15  gives  the  rebinding  of  Aq,  Aj  ,  and  Aj  at 
various  temperatures.  The  data,  taken  at  pH  6.8  in  75X  glycerol-water,  show 
the  following  features: 

(i)  At  all  temperatures,  the  four  substates  rebind  with  different  kinetics, 
(ii)  Rebinding  of  each  A  substate  is  nonexponential  in  time  up  to  180K. 

(ill)  A-j  is  always  considerably  slower  than  A[ . 

We  evaluate  the  nonexponential  kinetics  by  assuming  that  each  protein 
molecule  i «  in  a  substate  with  barrier  height  Hg^  for  rebinding  from  the 
pocket  state  B  to  the  covalently  bound  state  A,^  If  g(HgA)dHgA  denotes  the 
probability  of  finding  a  barrier  with  height  between  Hg^  and  Hg^  +  dHg^,  the 
rebinding  function  Is  given  by 
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N(t)  -  N(0)  /dHBA  g(HBA)  exp{-kBA(KBA)t}.  (3) 

Here  ''•ba^BA^  is  the  rate  coefficient  appropriate  for  a  barrier  of  height 
H0A.  Below  about  40K,  kBA^HBA^  18  dominated  by  tunneling.^’^  Above  40K, 
kBA(HBA^  *8  re^ate<^  t0  **y  an  Arrbenlus  relation, 

kBA^HBA^  "  ^BA^^O^  exPl-flBA',R^^  *  ^ 

We  have  verified  this  form,  with  the  overall  preexponential  factor  propor- 

28 

tional  to  T,  for  CO  binding  to  the  separated  6  chain  of  hemoglobin  Ziirich. 

Tq  is  an  arbitrary  reference  temperature  which  we  take  to  be  IOOK.  the 
parameter  AgA  consequently  gives  the  preexponential  factor  at  100K.  We 

further  assume  g(HgA^  to  1,6  8*ven  the  gamma  distribution  as  in  the  model  of 
29 

Young  and  Bowne.  This  distribution  peaks  at  the  activation  enthalpy 

Hpeak “  The  reblnding  parameters  are  found  by  fitting  the  experimental  data  to 

Eq.  (3),  with  g(HBA^  8lven  bY  the  gamma  distribution  and  kRA  by  Eq.  (4).  We 

5  28  30 

have  used  this  technique  for  many  heme  protein-ligand  systems.  ’  *  The 

barrier  parameters  for  binding  to  the  A  substates  are  given  in  Table  2.  For 
comparison,  we  also  list  the  parameters  obtained  by  monitoring  at  440  nm. 

3.4  The  B  Substates 

Since  we  observe  four  discrete  CO-bound  substates,  Aq-Aq,  we  also  expect 
four  discrete  stretching  frequencies  in  the  photodi ssoci ated  state  B.  We 
observe,  however,  only  three  bands,  as  shown  in  Fig.  lb.  In  an  earlier  paper 
we  have  discussed  these  bands  in  detail.  B ^  at  2119  cm  is  an  Intermediate 
that  decays  to  Bj  at  temperatures  above  15K.  For  the  treatment  of  rebinding 
above  about  15K,  we  are  left  with  only  Bq  and  Bp  At  pH  3.6  and  6.8,  B^  is 
much  more  Intense  than  Bq,  and  it  decays  with  the  same  time  dependence  as  Ap 
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while  Bq  shows  the  same  time  dependence  as  Aj.  B ^  therefore  must  bind  to 
substate  Aj  and  Bq  to  substate  A^. 

The  fact  that  we  cannot  see  the  B  substates  that  bind  to  Aq  and  Aj  1  s  not 

surprising.  The  Integrated  extinction  coefficients  of  the  B  bands  are  small 
—  1  -2 

(1.4  ±  0.1  mM  cm  )  and  a  sufficiently  sensitive  measurement  of  the  spectrum 
is  difficult.  Moreover,  the  intensities  of  Aq  and  A^  are  much  smaller  than 
that  of  Aj.  We  consequently  expect  the  B  substates  that  bind  to  Aq  and  A2  to 
be  much  smaller  than  Bj  and  therefore  to  be  hidden. 

3.5  Rebinding  Observed  Between  350  and  800  nm 

The  different  kinetics  for  the  various  A  substates  raises  the  question  of 
whether  similar  differences  can  also  be  seen  over  the  interval  350  to  800  nm, 
where  the  spectral  features  are  dominated  by  the  heme  group.  The  time  range 
from  100  ns  to  300  s  was  scanned  with  our  flash  photolysis  system  at  the 
wavelengths  indicated  in  Fig.  la.  Times  from  60  s  to  about  30  ks  were 
monitored  with  the  0 LIS-Cary  spectrometer  where  the  entire  wavelength  region 
was  scanned.  While  small  differences  were  observed,  no  major  variations  in 
the  time  course  as  a  function  of  wavelength  were  found  at  pH  6.8.  The  fastest 
rebinding  was  observed  around  the  maximum  of  the  CO  bound  spectrum,  the 
slowest  around  the  peak  of  deoxymyoglobin.  The  average  rates  at  421  and  440 
nm  differ  by  a  factor  of  2.2  at  60K  and  1.6  at  80K.  The  kinetics  at  440  nm, 
where  we  normally  monitor  rebinding,  agrees  well  with  that  of  the  integrated  A 
substates.  The  result  is  reassuring.  Observations  in  the  Soret  region  yield 
information  about  the  average  rebinding  behavior. 

3.6  Transitions  Among  A  Substates? 

We  have  shown  in  Subsection  3.1  that  the  A  substates  are  frozen  when  the 
protein  is  embedded  in  a  solid  matrix.  This  observation  leaves  open  the 
possibility  that  even  below  180K  transitions  among  B  substates  can  take  place. 
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Transitions  among  the  A  substates  could  then  be  Induced  by  light,  for  Instance 
through  the  chain  Aj  ♦  Bj  •*  B^  »  Aj.  To  study  this  possibility,  we  have 
performed  a  number  of  additional  experiments. 

(I)  Mb CO  After  Illumination.  If  transitions  among  the  B  substates  occur 
the  ratios  A^/Aj  measured  after  cooling  in  the  dark  and  after  illumination 
should  in  general  differ  below  180K.  We  have  performed  such  experiments  at 

1 60K,  140K,  and  120K.  The  change  in  relative  intensities  after  illumination 
is  less  than  1%  for  Aj  and  less  than  5Z  for  Aq  and  A^.  The  B  substates  do  not 
exchange  significantly  with  each  other  below  180K. 

(II)  Rebinding  Under  Illumination.  As  a  second  test,  we  measured  the  IR 
spectrum  as  a  function  of  time  during  exposure  of  the  MbCO  sample  to  constant 
illumination.  Under  the  assumption  that  each  A  substate  communicates  only 
with  one  B  substate  and  that  neither  A  nor  B  substates  exchange  with  each 
other,  the  rate  A  ♦  B,  denoted  by  k^,  is  determined  by  the  intensity  of 
illumination.  Rebinding  B  »  A  is  characterised  by  the  rate  k.  The  absorbance 
change  AAj(t)  in  the  band  Aj  is  then  given  by 

AAj(t)  “  AAj{kL/(kL+k)}{l-expl-(kL+k)t)},  (5) 

where  AAj  is  the  maximum  absorbance  change  of  band  Aj ,  corresponding  to 
complete  photodissociation.  If  the  rebinding  rate  k  is  given  by  a  distri¬ 
bution,  Eq.  (5)  must  be  generalized. 

For  kL  »  k,  Eq.  (5)  predicts  two  features:  (i)  The  approach  to  the 
steady  state  at  t  •  *  is  given  by  k^.  (ii)  The  steady-state  value  of  the 
absorbance  change,  AAj(“),  is  equal  to  AAj .  At  40K  both  of  these  predictions 
are  correct  to  within  about  3Z.  All  A  and  B  bands  approach  the  steady  state 


with  the  same  kinetics,  given  by  k^  and  the  steady  state  values  agree  with  the 
values  shown  in  Fig.  2.  This  result  verifies  that  no  transitions  Aj  -  A^  or 
Bj  -  Bj  take  place  at  AOK.  Data  at  60K  give  the  same  result.  At  70K  and 
above,  however,  Aj  approaches  equilibrium  much  slower  than  predicted  by  Eq. 
(5).  The  data  can  be  explained  if  a  small  fraction  of  the  Aj  substates  are 
pumped  into  longer-lived  states.  Observations  of  all  A  and  B  substates  rule 
out  transitions  from  Aj  to  other  A  substates.  We  discuss  an  explanation  of 
the  observed  effect  in  Subsection  A. 5. 

(iii)  Multiple-flash  Experiments.  As  a  third  test,  we  have  performed 
multiple-flash  experiments  at  various  temperatures.  In  these  experi¬ 
ments5,13’^,  photodissociation  is  repeated  before  all  Mb  molecules  have 
rebound  a  ligand.  Monitoring  of  the  separate  A  lines  can  reveal  if  pumping 
from  one  A  substate  to  another  occurs.  At  AOK,  flashes  were  repeated  every 
3000  s.  No  pumping  from  one  A  substate  to  another  is  observed  over  the  course 
of  A  flashes. 

(iv)  Photodissociation  After  Illumination.  In  the  earlier  multiple-flash 
experiments  (monitored  in  the  visible),  flashing  is  repeated  at  intervals  of 
length  r  on  the  order  of  20  times.3  If  transitions  from  one  substate  to  a 
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longer-lived  one  occurred  with  a  small  probability,  say  10  ,  or  only  at  times 

longer  than  T,  they  would  not  be  detected.  A  different  experiment  permits  the 
search  for  such  rare  transitions:  We  expose  the  sample  to  an  intense  white 
light  (kL  *  10  s-1)  for  a  time  tL  varying  between  10  s  and  103s.  The  light 
continuously  re-photolyzes  the  rebound  CO.  At  80K  and  above  more  than  50Z  of 
the  CO  molecules  rebind  within  less  than  0.1  s  and  will  therefore  be  re- 
photolyzed  between  10^  and  10^  times.  At  the  end  of  the  illumination,  the 
sample  is  completely  photolyzed  by  a  laser  pulse  and  rebinding  is  measured. 
Results  of  such  experiments  at  30K,  A0K,  92K  and  139K  are  given  in  Fig.  16. 
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At  30K,  a  very  small  Increase  In  long-lived  states  Is  seen  for  tL  -  1100  s. 

At  40K,  pumping  Is  still  small.  At  higher  temperatures,  rebinding  at  short 
times  remains  unaltered;  at  139K,  >90Z  of  the  photolyzed  CO  rebind  without 
change  In  rate.  The  change  In  the  rebinding  curve  at  long  times,  however.  Is 
dramatic;  long-lived  substates  are  greatly  enhanced. 

In  a  second  set  of  experiments,  a  second  laser  pulse  was  applied  with  a 
delay,  tD,  after  the  end  of  the  illumination.  Data  at  92K,  given  In  Fig.  17, 
prove  that  the  protein  relaxes  toward  the  initial  state  after  the  illumination 
is  terminated.  The  relaxation  Is  slow,  however,  and  Is  not  finished  at  time 
tjj  “  t^.  We  Interpret  these  results  in  terms  of  a  hierarchy  of  substates  in 


Subsection  4.5 
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4.  SUMMARY  AND  INTERPRETATION 

The  experiments  reported  In  Section  3  answer  the  main  questions  raised  in 
the  Introduction,  but  also  reveal  new  problems  and  point  to  new  avenues  of 
Investigation.  We  summarize  here  the  main  results  and  interpret  some  of  the 
new  data. 

4.1  Conformational  Substates 

One  of  the  motivations  for  the  present  work,  was  the  question  raised  by 
the  work  of  B.  Chance  et  al.^:  Is  rebinding  bi-  or  triphasic  as  originally 
claimed  by  Iizuka  et  al.^  or  is  it  characterized  by  a  power  law?'’  The  result 
of  the  present  work,  displayed  in  Figs.  14  and  15,  is  unambiguous:  Regardless 
of  where  rebinding  is  monitored,  the  kinetics  is  smooth  and  nonexponential 
over  more  than  12  orders  of  magnitude  in  time.  The  original  reason  for  intro¬ 
ducing  the  concept  of  conformational  substates  is  consequently  still  valid. 

The  observation  of  pumping  into  long-lived  states  (Figs.  16  and  17) 

raises,  however,  the  question  of  whether  the  arguments  rejecting  all  models 
32 

with  only  sequential  rather  than  parallel  barriers  are  still  valid.  The 
answer  is  yes.  The  multiple-flash  technique,  described  in  detail  else- 
where^’^’^  is  very  sensitive  to  the  difference  between  homogeneous  and 
inhomogeneous  systems  and  shows  that  the  main  components  of  the  system  must  be 
Inhomogeneous.  We  will  show  below  in  Subsection  4.5  (iv)  that  the  pumping 
experiments  lead  to  a  generalization  of  our  model  which  can  explain  all  data. 

4.2  The  Substates  in  the  Heme  Pocket 

(i)  The  A  Substates.  It  has  been  known  for  a  considerable  time  that  CO 

00  00  07  0/ 

binds  to  the  heme  iron  in  more  than  one  conformation.  ’  ’  ’  Indeed, 

Fig.  5  shows  four  different  stretching  bands  and  implies  four  substates,  Aq  to 
Aj.  The  reversible  exchange  among  the  substates  down  to  about  180K,  shown  in 
Fig.  6,  implies  that  all  bands  correspond  to  folded,  undenatured  proteins  and 
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can  occur  within  the  same  molecule.  The  main  properties  of  the  A  substates 
are  given  in  Table  1  and  Figs.  6-9. 

The  spectroscopic  data  in  Section  3  together  with  the  X-ray  structure 

35 

obtained  at  260K  by  Kuriyan  et  al.  permit  a  tentative  structural  assignment 
of  Aq  through  A^.  Kuriyan  et  al.  find  four  different  orientations  ("confor¬ 
mations  A,  B,  C,  D")  of  the  bound  CO  molecule  with  respect  to  the  heme  group, 
and  they  refine  two  of  the  four.  Their  conformation  C  has  an  occupation  of 
78%  and  can  therefore  be  identified  with  our  substate  A^  which  has  an 
occupation  of  70%  at  about  pH  6  and  260K.  With  this  identification,  Aj  has  an 

angle  Q(Fe-C-O)  of  141°  and  points  toward  C1D,  the  carbon  atom  1  of  pyrrole  D 
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of  the  heme  group.  (The  nomenclature  is  defined  in  Fig.  2  of  Kuriyan  et  al.  ) 
Their  conformation  D  has  an  occupation  of  22%.  Fig.  5  shows  that  Aj  has  a 
fractional  intensity  of  about  25%.  We  therefore  Identify  D  with  A-j.  A^  then 
has  an  angle  6  -  120°  and  points  in  the  direction  of  C4B,  a  considerable 
distance  away  from  Aj . 

We  identify  Aq  and  at  the  same  time  solve  an  old  problem  by  considering 
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the  effect  of  pH  on  Aq.  Watson  and  Kendrew  pointed  out  that  a  water 
molecule  bound  to  the  iron  in  metMb  would  be  within  hydrogen-bonding  distance 
of  the  distal  His  64(E7).  Since  CO  is  not  very  different  from  a  water 

molecule  in  size,  a  hydrogen  bond  between  the  oxygen  of  the  CO  and  His  E7 
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could  be  expected.  Neutron  diffraction  data  by  Hanson  and  Schoenborn  , 
however,  show  no  hydrogen  bond  in  MbCO.  If  Aq  is  identified  with  conformation 
A  of  Kuriyan  et  al.,  the  following  picture  emerges:  Of  the  four  substates, 

Aq  is  closest  to  His  E7  and  can  form  a  hydrogen  bond  with  Ne  if  His  E7  is 

protonated.  If  His  E7  is  the  heme-linked  protonation  group  with  pK  «  5.7 
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observed  in  optical  and  NMR  measurements  ’  we  expect  CO  in  the  substate  Aq 
to  be  fully  hydrogen-bonded  co  His  E7  well  below  pH  5.7  and  the  hydrogen  bond 
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to  be  absent  well  above  pH  5.7.  The  data  In  Table  1  and  In  Fig.  9  support 
this  expectation.  At  pH  5,  substate  Aq  Is  more  tightly  bound  than  Aj  by  about 
8  kJ/mol,  and  at  pH  6.8  It  Is  more  tightly  bound  by  about  A  kJ/mol.  The 
difference,  4  kJ/mol,  is  consistent  with  a  weak  hydrogen  bond.  Additional 
evidence  for  the  identification  of  Aq  comes  from  the  angle  0  *  154°  for 
conformation  A  which  is  similar  to  the  angle  found  In  chelated  protoheme- 
CO.^  Heme  carbonyls  typically^*  have  CO  stretching  frequencies  near  1970 
cm-1,  similar  to  the  value  for  Aq.  The  fact  that  no  hydrogen  bond  Is  seen  in 
the  neutron  diffraction  data  is  explained  by  the  low  Intensity  of  substate  Aq 
at  room  temperature  and  the  absence  of  protonated  Kls  E7  near  neutral  pH.  To 
test  the  assignments  directly,  the  Intensities  of  the  1R  bands  and  the  X-ray 
conformations  should  be  measured  on  the  same  crystal  at  a  pH  for  which  Aq  is 
sufficiently  populated.  Moreover,  additional  pH  studies  are  needed  to 
determine  whether  the  four  bands  correspond  to  four  distinct  orientations  or 
two  bands  correspond  to  the  protonated  and  deprotonated  forms  of  a  single 
orientation. 

We  summarize  the  complete  assignment  in  Table  3  where  we  also  include  the 
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van  der  Waals  energies  calculated  by  Kuriyan  et  al.  Table  3  shows  that  the 
calculated  binding  energies  are  close  to  the  relative  enthalpies  found  here 
but  are  generally  larger.  In  reality,  the  difference  in  occupation  is 
determined  as  much  by  the  relative  entropies  as  by  the  binding  enthalpies. 

(il)  Interpretation  of  the  B  Substates.  The  IR  spectrum  of  the  CO 
stretching  modes  in  the  photodlssoclated  state  B  is  shown  in  Fig.  lb.  Above 
20K  only  two  bands,  Bq  and  Bp  are  unambiguously  observed.  Bq  binds  to 
substate  and  Bj  to  Ap  At  present,  no  X-ray  diffraction  data  are  available 
on  the  photodlssoclated  CO,  so  no  sterlc  Interpretation  of  the  B  substates  can 
be  given.  The  fact  that  no  exchange  among  Bq  and  Bj  is  observed  below  about 
180K  suggests  that  each  is  sterlcally  close  to  the  corresponding  A  substate. 
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4.3  CO  Binding  at  Low  Temperatures 

(I)  Binding  Pathways.  Each  Mb  molecule  possesses  four  different  pathways 
within  the  pocket  for  CO  binding.  In  751  glycerol-water,  there  is  no  exchange 
among  the  four  channels  below  about  1 7 OK ;  above  180K  transitions  +  Aj  and 
presumably  Bj  ♦  B^  occur.  The  kinetics  of  each  of  the  four  pathways  is 
nonexponential  In  time  below  about  180K,  as  shown  in  Figs.  14  and  15.  Each  is 
governed  by  different  barrier  parameters  as  given  in  Table  2. 

The  data  provide  no  evidence  for  discrete  way-stations  or  "docking" 
within  the  heme  pocket  as  postulated  by  B.  Chance  et  al.4^  The  substate  B2 

Q 

decays  quickly  to  Bj  at  20K  and  above  ,  long  before  any  appreciable  rebinding 
occu rs •  B,  binds  to  substate  Aj  with  essentially  the  same  kinetics  as 
observed  in  the  visible.  Above  20K,  no  Intermediate  substates  appear  in 
either  the  bound  or  the  photodi ssoclated  CO  bands. 

(II)  Reaction  Theory.  The  data  in  Table  2  imply  that  both  entropic  and 
enthalpic  factors  control  the  binding  reaction  at  the  heme  iron.  A q  and  Aj 
both  have  about  the  same  peak  activation  enthalpy  H  a^;  the  faster  binding 
rate  of  Aq  relative  to  A^  is  caused  by  the  larger  preexponential  factor,  A^. 

Aj  has  a  slightly  larger  preexponential  than  Aj ,  but  has  a  slower  binding  rate 
owing  to  a  much  larger  peak  activation  enthalpy.  An  explanation  of  these 
characteristics  in  terms  of  the  pocket  structure  and  of  specific  residues  will 
require  detailed  studies  of  binding  rates  as  a  function  of  pH  and  with 
modified  myoglobins. 

The  parameters  listed  in  Tables  1  and  2  show  a  Brdnsted  correlation4^: 

The  most  tightly  bound  substate,  Aq,  binds  fastest  and  the  most  weakly  bound, 

A3,  binds  slowest.  The  same  correlation  is  found  for  different  heme  proteins 
where  the  one  with  the  most  tightly  bound  ligand  has  the  smallest  activation 
enthalpy  barrier.44  Such  a  correlation  is  expected  for  an  "overcoupled" 
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transition,  where  an  increase  in  the  well  depth  for  the  covalently  bound  state 
A  leads  to  a  decrease  in  the  barrier  between  the  pocket  state  B  and  A.  The 
experimental  result  consequently  demonstrates  that  CO  binding  to  many  heme 
proteins,  in  particular  Mb,  is  "overcoupled"  and  not  "undercoupled". ^ 

4.4  Influence  of  the  Environment  and  Control 

From  the  biological  point  of  view,  the  most  important  feature  of  our  data 

is  the  exquisite  sensitivity  of  the  occupation  of  the  substates  in  the  pocket 

to  external  parameters  such  as  pH  (Fig,  12),  pressure  (Fig.  11),  and  solvent 
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state  (Figs.  2-4).  ’  Earlier  we  found  that  the  association  rate  of  CO  to 

Mb  increases  with  decreasing  pH  and  follows  a  Henderson-Hasselbalch  equation 
with  pK  =■  5.7.  The  present  data  provide  more  insight  into  the  pH  dependence; 
a  change  in  pH  can  not  only  affect  the  kinetics  of  an  individual  pathway,  but 
can  also  shift  the  system  from  one  pathway  to  another.  Such  shifts  also  occur 
on  changes  in  pressure  and  solvent  state.  Three  conclusions  may  be  drawn  from 
these  experiments:  (i)  Processes  in  the  heme  pocket  are  extremely  responsive 
to  external  influer  es.  (ii)  IR  observation  of  ligand  binding  is  a  tool  well 
suited  for  the  investigation  of  the  effect  of  external  parameters.  (iii)  The 
sensitivity  of  the  processes  in  the  heme  pocket  to  external  parameters  may  be 
used  in  vivo  for  the  control  of  protein  reactions,  for  example  in  membrane 
proteins  through  changes  in  membrane  viscosity  and  composition,  and  in 
hemoglobin  through  quaternary  changes.  In  all  processes,  entropy  may  play  a 
more  important  role  than  enthalpy,  as  is  suggested  by  Table  2. 

4.5  Protein  Dynamics  and  Models 

The  sensitivity  of  the  IR  stretching  bands  to  the  conformation  of  the 
protein  makes  them  excellent  probes  for  investigating  protein  dynamics  and 
models.  We  sketch  the  background  for  such  investigations  and  describe  the 
conclusions  that  emerge  from  the  experiments  reported  in  Section  3. 
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(1)  The  Hierarchy  of  Conformational  Substates  (CS).  A  model  with  a 
single  set  of  CS  Is  not  capable  of  explaining  all  experiments.  On  the  one 

hand,  rebinding  in  a  75%  glycerol-water  solvent  is  nonexponential  up  to  at 
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least  170K  ,  and  pressure  titration  experiments  prove  that  transitions  among 

the  CS  responsible  for  the  nonexponential  kinetics  occur  only  above  about 
180K.^7  On  the  other  hand,  evidence  for  slow  transitions  among  substates  even 
at  80K  comes  from  X-ray  diffraction:  The  Debye-Waller  factor  of  many  atoms 
changes  appreciably  down  to  80K.^®*^  Studies  of  the  rearrangement  of  the 
protein  structure  after  photodissociation  ("proteinquake")  demonstrate  that  a 
hierarchy  of  motions  exists  and  that  protein  conformations  can  be  classified 
into  a  hierarchy  of  substates .  ^  We  draw  in  Fig.  18  the  hierarchy  of  CS  as 

modified  by  the  results  of  the  present  work. 

Figure  18  represents  the  conformational  energy  of  MbCO  drawn  for  various 
tiers  as  function  of  a  conformational  coordinate  (cc).  The  energy  valley  in 
the  top  diagram  represents  MbCO.  The  observation  of  four  CO  stretching  bands 
implies  that  MbCO  can  exist  in  four  distinct  substates  of  the  zeroth  tier 
(CS^),  Aq  to  Aj.  The  nonexponential  rebinding  kinetics  observed  for  each  band 
suggests  that  the  actual  energy  surface  for  each  CS^*  consists  of  a  large 
number  of  conf ormational  substates  of  the  first  tier,  CS^.  Each  of  these 
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again  is  subdivided  into  a  large  number  of  substates  of  the  second  tier,  CS  . 
The  furcation  continues,  but  the  lower  tiers  are  not  relevant  for  the  present 
discussion. 

At  rest,  a  protein  will  not  remain  in  one  CS,  but  equilibrium 
fluctuations  (EF)  will  move  the  protein  from  CS  to  CS.  The  rates  of  the  EF 
depend  strongly  on  temperature  and  environment.  In  75%  glycerol-water,  only 
EF2  (and  EF3...)  take  place  below  about  170K,  but  each  protein  remains  frozen 
in  a  particular  CS1.  Above  about  180K,  EF1  occur  and  the  protein  fluctuates 


rapidly  from  one  CS^  to  another.  After  photodissociation  of  MbCO,  proteins 
relax  from  one  state  (MbCO)  to  another  (Mb)  through  a  sequence  of  functionally 
Important  motions  (fims).^ 

(11)  Relaxation  of  the  A  Substates.  In  a  75Z  glycerol-water  solvent,  the 
behavior  of  the  A  substates  undergoes  a  change  near  180K.  The  CS®  no  longer 
Interconvert  (Fig.  6),  the  band  frequency  Vqq  shifts,  and  the  linewidth 
changes.  We  interpret  these  features  in  terms  of  the  CS1.  If  the  large-scale 
motions  of  the  proteins  are  absent,  the  controlling  stdechains  within  the 
pocket  may  be  prevented  from  fluctuating,  thereby  stopping  transitions  among 
the  CS^.  This  explanation  also  accounts  for  the  large  change  in  the  linewidth 
T  of  band  Aq  in  Fig.  9.  At  pH  5  almost  all  of  the  His  E7  are  protonated  and 
therefore  most  of  the  CO  in  state  Aq  are  hydrogen  bonded.  Below  about  180K, 
the  motion  of  His  E7  and  therefore  also  the  motion  of  the  CO  is  restricted, 
leading  to  a  narrow  band.  Above  18QK,  the  sidechain  motion  may  also  force  the 
CO  to  large  excursions,  leading  to  a  large  T.  At  pH  6.8,  approximately  10-20Z 
of  the  CO  are  hydrogen  bonded  to  His  E7  and  these  may  account  for  the  smaller 
increase  in  r  above  180K. 

(ill)  Slaved  Glass  Transition.  Mb  in  75%  glycerol-water  possesses  a 
remarkable  transition  between  about  180  and  200K  as  discussed  in  (il).  The 
transition  occurs  in  an  interval  of  about  20K  (Figs.  6  and  8;  and  approach  to 
equilibrium  is  nonexponential  in  time  (Fig.  10).  These  characteristics  are 
similar  to  the  ones  observed  in  glass  and  spin  glass  transitions.  Glasses 
(amorphous  solids )"* ^  ^ ^  and  spin  glasses'*4  are  disordered  and  frustrated'’"’ 
systems  that  are  believed  to  have  a  highly  degenerate  (many-valley )  ground 

state. ^  The  transition  to  the  glass  or  spin-glass  state  occurs  in  a  finite 
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temperature  range  and  Is  nonexponential  in  time.  ’  Because  proteins  are 

18 

disordered  and  most  likely  frustrated  systems  with  many  energy  valleys 
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(Fig.  18)  we  expect  a  transition  similar  to  a  glass  or  spin-glass  transition. 
The  transition  near  180K  satisfies  all  the  requl resents . 

This  transition  manifests  one  property  not  encountered  in  glarses:  It 
depends  crucially  on  the  solvent  surrounding  the  protein.  In  75Z  glycerol- 
water,  the  transition  is  at  approximately  180K  (Fig.  8).  The  solvent  undergoes 
a  glass  transition  at  essentially  the  same  temperature  (Fig.  13). 5^’^®  In 
water,  the  transition  Is  near  260K,  close  to  the  melting  point  of  .he  buffered 
solvent  (Fig.  6).  In  solid  PVA,  no  transition  in  MbCO  is  apparent  up  to  at 
least  300K,  although  the  ratio  Aj/A^  decreases  gradually  above  100K.  The 
observations  Imply  that  the  EF1,  the  fluctuations  of  the  CS^,  are  suppressed 
when  the  protein  surrounding  is  frozen  or  solid.  The  motions  of  the  first 
tier  thus  are  "slaved"  to  the  motions  of  the  solvent  and  we  consequently  call 
the  transition  a  "slaved  glass  transition." 

The  slaved  glass  transition  is  not  a  property  of  the  protein  alone; 
protein  and  solvent  together  must  be  considered  as  one  system.  The  hydration 
shell  of  the  protein  most  likely  plays  a  crucial  role  both  in  the  glass 
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transition  and  in  the  function  of  the  protein.  *  ’  Studies  of  the 
hydration  water  in  Mb  solutions  and  crystals  by  calorimetry  and  IR  spectros¬ 
copy  show  a  broad  glass  transition  between  180  and  270K  which  depends  on  the 
degree  of  hydration.^5 

(iv)  Pumping  and  Substate  Symmetry  Breaking.  The  results  in  Figs.  16  and 
17  appear  at  first  to  call  for  a  new  explanation  of  the  nonexponential 
kinetics.  We  show  here,  however,  that  an  extension  of  our  earlier  model 
explains  the  pumping  into  long-lived  states.  In  that  model5  we  ascribe  the 
nonexponential  kinetics  to  a  distribution  of  activation  enthalpies,  charac¬ 
terized  by  the  probability  density  g(HBA)  as  given  for  the  binding  of  CO  to  Mb 
in  Fig.  19a.  Each  substate  CS^  of  the  first  tier  in  Fig.  18  is  characterized 
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by  a  unique  HgA»  Different  CS  ^  thus  have  different  properties:  the  energy 
valleys  of  the  first  tier  are  not  identical  and  we  call  this  property 
"substate  symmetry  breaking".  At  this  level  of  substate  symmetry  breaking  the 

2 

substates  of  the  second  tier,  CS  ,  have  no  effect  on  the  activation  energy  for 
covalent  bond  formation.  Experiments  suggest  that  this  assumption  may  be 
incorrect.  The  nonexponential  relaxation  in  tier  2  seen  in  the  protelnquake 
experiment^  and  the  temperature  dependence  of  the  Debye-Waller  factor  below 
180K^®’^  imply  that  different  CS^  have  different  physical  properties.  The 

substate  symmetry  thus  is  also  broken  in  the  second  tier,  and  we  assume  that 

2  1 
different  CS  have  different  barriers  for  the  binding  step  B  +  A.  The  CS 

2 

determine  the  barrier  crudely,  while  the  CS  fine-tune  the  barrier.  We 
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represent  the  situation  in  Fig.  19b.  Within  each  CS  ,  the  CS  determine  a 
range  of  activation  enthalpies.  The  distribution  of  the  barriers  within  a 
given  CS1  is  indicated  by  the  solid  curves,  and  the  overall  distribution  as 
determined  by  the  CS1  is  given  by  the  dashed  envelope. 

A  possible  explanation  of  the  pumping  observed  in  Figs.  16  and  17  is  now 
straightforward.  Since  the  CS^  are  frozen  below  180K,  we  only  have  to 
consider  a  given  CS^.  Fig.  20  gives  a  schematic  representation  of  the 
barriers  for  the  binding  step  B  ♦  A  as  a  function  of  a  conformational 
coordinate  cc2  of  the  second  tier.  The  two  coordinates  have  a  very  different 
meaning.  The  conformational  coordinate  cc2  describes  the  protein  structure 
within  a  given  CS^;  shown  are  three  substates  of  the  second  tier,  labelled  by 
a,  b,  and  c.  The  reaction  coordinate  rc  describes  the  transition  of  the  CO 
from  the  pocket  (B)  to  the  covalently  bound  state  at  the  iron  (A).  The 

potentials  along  the  reaction  coordinate  for  fixed  cc2  are  drawn  as  dashed 

2 

lines.  If  the  system  is  initially  in  a,  the  lowest  CS  ,  photodissociation 
will  move  it  to  a'.  The  system  can  now  either  rebind,  B  ♦  A,  without  change 
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2 

in  cc2,  or  it  can  relax  to  another  CS  of  lower  Gibbs  energy,  for  instance  b'. 

2 

We  have  studied  the  relaxation  to  lower  CS  through  flm  2  in  detail  and  have 
found  that  it  is  nonexponential  in  time  and  extends  over  extremely  large 
ranges  in  time.17  In  states  b'  and  c' ,  the  barrier  B  ♦  A  may  either  be  higher 
or  lower  than  in  a'.  If  it  is  lower,  the  system  will  rapidly  return  to  b  and 
from  there  relax  to  a.  If  it  is  higher,  the  system  will  stay  longer  in  B  and 

this  behavior  appears  as  pumping.  The  shift  to  longer  rebinding  times  is 

64 

caused  by  a  relaxation  in  the  protein  coordinate  and  not  by  a  longer-lived 
Intermediate  state  between  B  and  A.  The  model  explains  both  the  pumping  seen 
in  Fig.  16  and  the  relaxation  in  Fig.  17.  Relaxation  occurs  if  the  system  is 

2 

first  pumped  into  long-lived  CS  ,  some  of  which  rebind  in  the  dark,  from,  say, 

2 

c'  to  c.  From  there  it  will  relax  toward  a,  the  lowest  CS  .  A  flash  at  a 
later  time  explores  the  degree  of  relaxation. 

A  theory  of  pumping  remains  to  be  constructed,  but  simple  estimates 
indicate  that  the  relaxation  times  and  enthalpies  observed  are  consistent  with 
our  earlier  experiments.17  Consider  the  pumping  observed  at  30K  in  Fig.  16. 
The  absorbance  scale  is  greatly  expanded.  The  rate  at  N(t)  «  0,95  changes  by 
a  factor  of  2  after  illuminating  for  1100  s.  This  rate  change  corresponds  to 
a  shift  of  6H®  *  RT  In  2  *  0.17  kJ/mol.  Since  8(Hg^)  extends  from  5  to  25 
kJ/mol,  such  a  small  shift  is  not  even  visible  in  Fig.  19  and  it  is  reasonable 
to  attribute  it  to  substates  of  the  second  tier.  At  139K,  the  largest  time- 
shift  in  Fig.  16  corresponds  to  an  enthalpy  shift  of  SHg^  »  8  kJ/mol.  This 
shift  will  change  the  distribution,  but  is  still  reasonable.  Times  and 
temperatures  seen  in  pumping  experiments  are  also  consistent  with  flm  2 
values.  Recent  measurements  on  fim  2  (unpublished)  indicate  that  it  starts 
already  at  30K  and  extends  to  ouch  longer  times  than  shown  in  ref.  17. 
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4.6  Remarks  and  Outlook 

(1)  The  IR  Bands  as  Tools.  The  experiments  described  and  evaluated  here 
demonstrate  that  IR  spectra,  when  studied  as  a  function  of  time,  temperature, 
solvent,  and  pressure,  provide  a  powerful  tool  for  studies  of  protein  dynamics 
and  protein  function.  Nearly  all  investigations  here  mark  a  beginning,  and 
extensions  to  shorter  times,  other  systems,  and  other  solvents  are  likely  to 
yield  even  richer  information. 

(il)  The  Chance  Experiment.  B.  Chance  et  al.^  have  made  two  claims, 
that  CO  rebinding  to  Mb  is  bi-  or  t ri exponential  at  40K  and  that  they  can  pump 
short-lived  into  long-lived  states.  Fig.  14  proves  that  the  first  claim  is 
incorrect;  data  extended  to  10^s  show  no  evidence  for  Individual  components. 
The  second  claim  has  to  be  discussed  in  more  detail  because  of  our  new 
evidence  for  a  small  amount  of  pumping  at  40K  in  Fig.  16.  Chance  et  al.  use  a 
very  weak  pump  light,  adjusted  so  that  on  the  average  about  50%  of  the  MbCO 
are  photolyzed.  Fig.  14  shows  that  N(t)  at  40K  drops  50%  by  about  5*10^s. 

Thus  in  the  Chance  experiment,  equilibrium  can  occur  only  for  illumination 

4 

times  much  greater  than  5«10  s.  Owing  to  the  large  distribution  of  flashoff 
rates  k^  through  the  thick  sample,  a  time-dependent,  nonequilibrium 
distribution  of  recombination  rates  is  established  in  the  sample.  In  the 
Chance  et  al.  experiment,  the  relative  amplitude  of  the  slow  components  is 
enhanced  by  up  to  a  factor  of  two  resulting  in  apparent  pumping.  Their 
result,  therefore,  can  be  explained  without  pumping  in  a  straightforward  way; 
it  is  caused  by  the  combination  of  weak  pump  light,  thick  sample, 
nonexponential  kinetics,  and  adjustment  to  about  50%  photolysis.  A  computer 
simulation  of  this  combination,  based  on  the  generalization  of  Eq.  (5), 
reproduces  the  apparent  pumping. 
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(Hi)  Theories.  A  goal  of  protein  dynamics  studies  is  a  theory  that 
classifies  and  describes  protein  states  and  protein  motions.  Remarkable 
progress  toward  this  goal  has  been  made  with  molecular  dynamics  computations , 
particularly  by  Karplus  and  collaborators.6^  Elber  and  Karplus66  have 
verified  the  existence  of  multiple  energy  valleys  with  a  300  ps  simulation. 

One  unsolved  question  in  such  calculations  concerns  the  level  of  substates 
that  are  reached.  Comparison  of  the  timescales  simulated  with  the  ones 
involved  in  the  various  tiers  of  substates17  suggests  that  molecular  dynamics 
explores  many  of  the  substates  of  tier  2,  but  cannot  yet  reach  tier  1. 

A  second  theoretical  approach  is  based  on  ultrametricity."16  We  have 
proposed  that  the  arrangement  of  the  substates  as  shown  in  Fig.  18  suggests 
ultrametricity. 17  Elber  and  Karplus  have  approached  this  problem  with 
molecular  dynamics66  and  find  no  evidence  for  ultrametricity.  Since  they 
investigate  the  distance  matrix  and  we  propose  ultrametricity  in  the  space  of 
conformational  energy,  no  disagreement  between  the  two  models  exists  yet.  One 
possible  problem  regarding  ultrametricity  is  raised  by  Fig.  18:  We  have 
assigned  the  zeroth  tier,  CS®,  to  the  four  A  substates.  Such  an  assignment  is 
certainly  correct  below  the  glass  temperature  where  no  exchanges  among  the  A 
substates  take  place.  Above  the  glass  temperature,  however,  both  the  CS11  and 
CS1  interchange,  and  it  is  no  longer  clear  how  to  classify  and  order  these 
substates.  Additional  investigations,  both  experimental  and  theoretical,  may 
shed  new  light  on  these  aspects  of  protein  states  and  motions. 
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TABLE  1 

Relative  binding  parameters  of  the  A  substates. 

R  Is  the  gas  constant.  Solvent:  75Z  glycerol/water. 


Substate 


*0 

A1 

*3 


HA(kJ/mol> 
pH  6.8  5. 

-A  -8 

0  0 


2 


sa/r 


6.8  5. 

-5  -7 

0  0 


2 


0 


0 
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TABLE  2 


Barrier  parameters  for  CO  binding  to  A  Substates. 


Substate 

Hpeak 

kJ/mol 

logIABA(100K)/s_l 

*0 

10 

10.8 

A1 

9.5 

9.3 

a3 

18 

9.8 

Soret  (440nm) 

10.1 

9.0 
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TABLE  3 


The  values  for  vCq  refer  to  a  752  glycerol-water  solvent,  pH  6.8, 

and  240K.  The  conformations  and  their  properties  are  from 
35 

Kuriyan  et  al.  ,  at  260K  and  pH  6.  4  and  0  characterize  the 

CO  conformation;  *  is  the  dihedral  angle  between  the  Fe-C-0 
and  the  Fe-C-NC  planes;  0  is  the  angle  Fe-C-O. 33 


Substate 

VC0 

ha 

Conformation 

F  35 

vdW 

4 

0 

-1 

cm 

kJ/mol 

kJ/mol 

1966 

-u 

A 

-12 

20° 

154 

Al 

1945 

0 

C 

-12 

60° 

141 

a2 

1941 

- 

(B) 

+25 

93“ 

120' 

A3 

1930 

2 

D 

+4 

O 

CM 

O 

1 

120 
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FIGURE  CAPTIONS 

1  The  spectral  range  where  rebinding  is  monitored.  (a)  shows  the 
absorption  spectra  of  Mb  (dotted  line)  and  MbCO  (solid  line)  in  the 
Soret,  the  a  and  8,  and  the  near-IR  regions.  Arrows  indicate  where  we 
have  measured  rebinding  in  the  flash  photolysis  experiments,  (b)  gives 
the  IR  difference  spectrum  (Mb-MbCO). 

2  Temperature-dependent  IR  spectrum  of  the  CO  stretching  bands  of  MbCO. 
Solvent:  755!  glycerol-water;  pH  6.8. 

3  Temperature-dependent  IR  spectrum  of  the  CO  stretching  bands  in  MbCO. 
Solvent:  water;  pH  6.7. 

4  Temperature-dependent  IR  spectrum  of  the  CO  stretching  bands  in  MbCO. 
Solvent:  solid  poly(vinyl)alcohol  (PVA);  pH  7.0. 

5  Voigtian  decomposition  of  the  MbCO  IR  spectrum  into  four  A  substates  at 
10K  and  240K.  Solvent:  75Z  glycerol-water;  pH  6.8. 

6  Plots  of  (a)  (Aq/Aj)  and  (b)  A^/Aj  versus  10^/T.  Aj/Aj  is  the  ratio  of 
the  areas  of  the  stretching  bands  in  the  bound  states  Aj  and  Aj.  Below 
180K,  the  ratios  change  little. 

7  Temperature-dependent  area  of  the  CO  stretching  bands  in  MbCO. 

8  Temperature  dependence  of  the  peak  frequency  of  the  CO  stretching  bands 
Aq  and  Aj  in  755!  glycerol-water  at  pH  5.0  and  pH  6.8. 
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9  Temperature  dependence  of  the  line  widths  (full  width  at  half  maximum)  of 
the  A  substates.  Solvent:  75Z  glycerol-water;  pH  6.8. 

10  Approach  to  equilibrium  of  substate  Aq  after  a  change  in  temperature  from 
195K  to  190K.  ng(t)  denotes  the  fraction  of  MbCO  that  have  not  reached 
equilibrium  in  substate  Aq  at  the  time  t  after  the  temperature  change. 

11  The  CO  stretching  frequencies  of  MbCO  at  0.1  MPa,  50  MPa,  100  MPa 

(1  MPa  »  9.9  atm).  Solvent:  75X  glycerol-water;  0.25  M  Tris-HCl  buffer; 

pH  6.1. 

12  The  CO  stretching  frequencies  of  MbCO  in  75Z  glycerol-water  at  various 
values  of  pH. 

13  1R  spectrum  of  75Z  glycerol-water  solvent  as  a  function  of  temperature. 

14  Rebinding  Mb  +  CO  ♦  MbCO  at  40K.  O  :  440  nm;  9:  Integrated  areas  of 
the  A  bands;  0:  Integrated  area  in  the  Soret  region,  400-450  nm; 

O  :  Integrated  area  of  the  760  nm  band.  The  individual  rebinding  curves 
have  not  been  matched  but  have  been  normalized  independently. 

15  Rebinding  to  the  substate  (a)  Aq,  (b)  Ap  and  (c)  A-j  as  a  function  of 
time  at  various  temperatures.  •  :  mld-lR  flash  photolysis  data. 

H  :  FTIR  data.  The  solid  lines  are  the  fits  to  the  data  using  Eqs.  (3) 
and  (4)  and  a  gamma  distribution  for  g(Hg*)*  The  fit  parameters  Hpeak 
and  AgA  are  listed  in  Table  2.  Solvent:  75X  glycerol-water;  pH  6.8. 
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16  Rebinding  after  Illumination.  The  MbCO  sample  is  Illuminated  by  an 
intense  white  light  for  tL  seconds  before  being  photodissoclated  by  a 
laser  flash.  Note  the  expanded  scale  for  the  absorbance  change  in  (a). 

In  (b)  the  lines  have  been  drawn  through  each  data  point.  Marker  symbols 
have  been  omitted  to  avoid  clutter. 

17  Photodissociation  after  a  delay  t q  following  illumination  for  t^  *  1110s. 
Normal  92K  rebinding  kinetics  (t^  “0,  tp  *  0)  is  included  for  comparison. 
Lines  have  been  drawn  through  each  data  point.  Solvent:  752  glycerol- 
water;  pH  7.0. 

18  Hierarchical  arrangement  of  conformational  substates  (CS)  in  Mb.  The 

conformational  energy  G  is  drawn  as  a  function  of  a  conformational 

3  4 

coordinate.  The  two  lowest  tiers,  CS  and  CS  ,  are  omitted. 

19  (a)  Activation  enthalpy  distribution  g(HgA)  determined  entirely  by 
conformational  substates  of  the  first  tier,  (b)  The  overall  distribu¬ 
tion,  the  dashed  envelope,  is  determined  by  the  substates  of  the  first 
tier;  the  substates  of  the  second  tier  provide  the  fine-tuning  represented 
by  the  narrow  solid  distributions. 

20  Rebinding  within  one  conformational  substate  of  the  first  tier.  The 
diagram  shows  the  energy  surfaces  as  a  function  of  the  conformational 
coordinate  cc2  and  the  reaction  coordinate  rc.  Details  are  given  in  the 
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